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Coverage of space in Boolean models 
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Abstract: For a marked point process {(xi, Si)i>\} with {xi £ A : i > 1} 
being a point process on A C M. d and {Si C R d : i > 1} being random sets 
consider the region C = Ui>\(xi + Si). This is the covered region obtained 
from the Boolean model {(xi + Si) : i > 1}. The Boolean model is said to 
be completely covered if A C C almost surely. If A is an infinite set such that 
s + A C A for all s £ A (e.g. the orthant), then the Boolean model is said to 
be eventually covered if t + A C C for some t almost surely. We discuss the 
issues of coverage when A is R d and when A is [0, oo) d . 



1. Introduction 

A question of interest in geometric probability and stochastic geometry is that of 
the complete coverage of a given region by smaller random sets. This study was 
initiated in the late 1950's. An account of the work done during that period may be 
found in Kendall and Moran (1963). A similar question is that of the connectedness 
of a random graph when two vertices u and v are connected with a probability p u — v 
independent of other pairs of vertices. Grimmett, Keane and Marstrand (1984) and 
Kalikow and Weiss (1988) have shown that barring the 'periodic' cases, the graph 
is almost surely connected if and only if = oo. 

Mandelbrot (1972) introduced the terminology interval processes to study ques- 
tions of coverage of the real line R by random intervals, and Shepp (1972) showed 
that if S is an inhomogeneous Poisson point process on R x [0, oo) with density 
measure A x fj, where A is the Lebesgue measure on the ir-axis and fi is a given 
measure on the y-axis, then \J^ X y\^g(x,x + y) = R almost surely if and only if 
J Q dx exp(J^°(y — x)ii(dy)) = oo. Shepp also considered random Cantor sets de- 
fined as follows: let 1 > t\ > ti > . . . be a sequence of positive numbers decreasing 
to and let V\,V%,... be Poisson point processes on R, each with density A. The 
set V := R \ (Uj U X £Vi {x,x + ti)) is the random Cantor set. He showed that V 
has Lebesgue measure if and only if J^i^i = 00 ■ Moreover, P(V = 0) = or 1 
according as J2n°=i n ~ 2 exp{A(<i + • • • + t n )} converges or diverges. 

In recent years the study has been re-initiated in light of its connection to percola- 
tion theory. Here we have a marked point process {ixi, Si)i>\} with {xi : i > 1} be- 
ing a point process on A C R d and Si C R d being random sets. Let C = Ui>i(xi+Si) 
be the covered region of the Boolean model {(xi + Si) : i > 1}. 

The simplest model to consider is the Poisson Boolean model, i.e., the process 
{xi : i > 1} is a stationary Poisson point process of intensity A on R d and 

^ = [0, Pl ] d , »>1 (1) 



*This research is supported in part by a grant from DST. 

Statistics and Mathematics Unit Indian Statistical Institute, 7 SJS Sansanwal Marg, New 
Delhi 110016, India, e-mail: rahulOisid. ac . in 

AMS 2000 subject classifications: primary 05C80, 05C40; secondary 60K35. 
Keywords and phrases: Poisson process, Boolean model, coverage. 



119 



120 



R. Roy 



are d-dimensional cubes the lengths of whose sides form an independent i.i.d collec- 
tion {pi : i > 1} of positive random variables. (Alternately, Si's are <i-dimensional 
spheres with random radius pi.) 

In this case, Hall (1988) showed that 

Theorem 1.1. C = M. d almost surely if and only if Epf = oo. 

More generally, Meester and Roy (1996) obtained 

Theorem 1.2. For {xi : i > 1} a stationary point process, if Epf = oo then 
C = R. d almost surely. 

The above results relate to the question of complete coverage of the space M. d . 

Another question which arises naturally in the Poisson Boolean model is that 
of eventual coverage (see Athreya, Roy and Sarkar [2004]). Let {xi : i > 1} be 
a stationary Poisson process of intensity A on the orthant and the Boolean 
model is constructed with random squares Si as above yielding the covered region 
C = Ui>i[xi(l),a:<(l) +pi] x ■ • • x [xi(d), Xi{d) + pi]. In that case, P(R^ C C) = 0, 
however we may say that is eventually covered if there exists < t < oo such 
that (t, oo) d C C. In this case, there is a dichotomy vis-a-vis dimensions in the 
coverage properties. In particular, while eventual coverage depends on the intensity 
A for d = 1, for d > 2 there is no such dependence. 

Theorem 1.3. For d = 1, 

(a) if < I := liminf^—xx, xP(p± > x) < oo then there exists < Ao < j < oo 
such that 



P\(R+ is eventually covered by C) 



if A < A 

1 if A > A ; 



(b) if < L := limsup^.^^ xP{p\ > x) < oo then there exists < j- < Ai < oo 
such that 



P\(R+ is eventually covered by C) 



if A < Ai 
if A> Ai; 



(c) if linx^oo xP(pi > x) ~ oo then for all A > 0, R+ is eventually covered by 
C almost surely (P\); 

(d) z/lim^^oo xP(pi > x) — then for any A > 0, R+ is not eventually covered 
by C almost surely (P\). 



Theorem 1.4. Let d>2, for all A > 
if 



(a) Pa(R+ is eventually covered by C) — 1 whenever\immf x ^ OQ xP(pi > x) > 



0. 

(b) Pa(K+ is eventually covered by C) — whenever lim^^oo xP{p\ > x) = 0. 

In 1-dimension for the discrete case we may consider a Markov model as follows: 
X\,X<2,... is a {0, 1} valued Markov chain and Si := [0, /?»], i = 1,2,... are i.i.d. 
intervals where is as employed in Q). The region U(i + Si)lx i= i is the covered 
region. This model has an interesting application in genomic sequencing (Ewens and 
Grant [2001]). If pij ~ P(X n+ i = j | X n = i), i,j = or 1, denote the transition 
probabilities of the Markov chain then we have 

Theorem 1.5. Suppose < Poo,Pw < 1- 
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(a) If I := liminfj^oo jP(pi > j) > 1, then P{C eventually covers N} = 1 when- 
ever P'! 1 > III. 

PlO+POl ' 

(b) If L := limsupj^,^, jP(fti > j) < oo, then P{C eventually covers N} = 
whenever — 291 — < 1/L. 

Molchanov and Scherbakov (2003) considered the case when the Boolean model 
is non- stationary . For a Poisson point process {xi : i > 1} of intensity A, we place 
a (i-dimensional ball B(xt, pih(\\xi\\)) centred at Xi and of radius /}j/i(||xj||) where 
pi is as before and h : [0, oo) — > (0, oo) is a nondecreasing function. Let 

C = U£ 1 B(x itPi h(\\x i \\)) 

denote the covered region. Let TTd denote the volume of a ball of unit radius in d 
dimensions and take ho (r) — ( j^- log r J 

Theorem 1.6. Suppose Ep d+r> < oo for some r\ > and h is as above. 
(a)// 

l h := liminf ( > — then P(C = R d ) > 0, 



h (r)J E(pf) 



h(r)\ d 1 



and 
(b) if 



0<£ft:=limsup hrTT <771i then P ( C = R ) = °- 

The result in (a) above cannot be translated into an almost sure result because 
of the lack of ergodicity in the model. 



2. Complete coverage 

We now sketch the proofs of Theorems 11.11 11.21 and 11.61 

Let V = [0, l] d \ C denote the 'vacant' region in the unit cube [0, l] d ; 

E(£(V)) = E l {x ig not covere d} dx 

J[o,i] d 

= exp(-XE(4) (2) 

where I stands for the d-dimensional Lebesgue measure. Hence, if Epf < oo then 
E(£(V)) > and so P([0, l} d C C) < 1. 

Conversely, Epf = oo implies E(£(V)) = and thus by stationarity, E(£(R d \ 
C) ) =0. Using the convexity of the shapes Si we may conclude that P(C = M. d ) = 1. 

Here the Poisson structure was used to obtain the expression (J2J; for a general 
process we need to extract, if possible, an ergodic component of the process and 
show that the Boolean model obtained from this ergodic component covers the 
entire space when Epf — oo. To this end let {xi : i > 1} be an ergodic point 
process with density 1. Let D n = [0, 2 n l d \ d and E n = {there exists xi in the 
annulus D n+ i \ D n such that Dq C (x, + Si)}. Also let A m be the event that m is 
the first index such that #{z : x, € -D ?i +i \ D n } > a2 n for all n > m and for some 
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fixed constant a. By ergodicity, {A m : m > 1} forms a partition of the probability 
space and we obtain 

P(n%L m E c k I A m ) < P(n? =m n i:XieDk+lXDk { Pl < 2 k ' d + 1} | A rn ) 

< n ^ 2k/d + i ) a2k+1 



k—m 

oo 



k— m+1 



< n n p(pf<^ k +j) 

k=m+l \ j=l 



n (i-p( P f>k)) 

k=2 m 

oo 

= if and only if P(pf > A:) = oo. 

fc=2 m 

This completes the proof of Theorems 11.11 and 11.21 

To prove Theorem ll.fcil fa) we study the case when h(r) = l^/ d ho{r) and IhEpf > 
1 + S for some 5 > 0. It may be easily seen that for z e Z d , P{z + (-1/2, l/2] d g 
C} < cxp{-A^(i? 2 )}, where i? z = {(r,a;) € [0, oo) x M d : z + (-1/2, l/2] d C 
B(x, rh{x)) and /it is the product measure of the measure governing pi and Lebesgue 
measure. From the properties of ho it may be seen, after some calculations, that 
given e > 0, there exists r such that for ||z|| > r, p(R z ) > (1 — e)i^d{h{z)) d Epf. 
Thus we obtain, for some constant K, 

£ P{z + (-l/2,l/2] d £C} < Kj2^M-d(l + 5)(l-e)log(\\z\\)} 



zll - d (l +S){ l-e) <00 



whenever (1 + — e) > 1. Invoking the Borel-Cantclli lemma we have that 
z + (— l/2,l/2] d % C occurs for only finitely many z G Z d . Using this we now 
complete the proof of Theorem 1 1.61 (a). 

The proof of Theorem ll.6r b) is more delicate and we just present the idea here. 
For d > 2, if we place points in a spherical shell of radius n 7 , such that the interpoint 
distances are maximum and are of the order of nP where < (3 < 7 and 7 > 1 then 
the number of points one can place on this shell is of the order of n 7 ^'^ 1 ' . Let V n 
be the event that one such point is not covered by C. It may be shown that there 
is a choice of 7 and (3 such that infinitely many events V n occur with probability 1. 
For d = 1, the same idea may be used and, in fact, the proof is much simpler. 



3. Eventual coverage 

We discretise the space by partitioning it into unit cells {{ii, ■ ■ . ,id) + (0, l] d : 
ii, . . . , id — 0, 1, . . .} and call a vertex i := (ii, . . . , id) green if x G i + (0, l] d for 
some point x of the point process. Consider two independent i.i.d. collections of 
random variables {p^} and {p\} where the distribution of p\ and p\ are identical 
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to that of 2 + Lmaxjpi, . . . , pn}\ and max{0, |_max{pi, . . . , pn}\ — 1} respectively; 
here N is an independent Poisson random variable with mean A conditioned to be 
1 or more. We now define two discrete models, an upper model and a lower model: 
in both these models a vertex i is open or closed independently of other vertices, 
and the covered region for the upper model is U{j open}(i + [0,/^]), and that for 
the lower model is U{i open}(i + [O,^])- Observe that the eventual coverage of the 
Poisson model ensures the same for the upper model and eventual coverage of the 
lower model ensures the same for the Poisson model. Thus it suffices to consider the 
eventual coverage question for a discrete model, as in Proposition 13.11 below, and 
check that the random variables p u and p l satisfy the conditions of the proposition. 

We take \X\ : i € N d } to be an i.i.d. collection of {0, 1} valued random variables 
with p := P(X\ — 1) and {p\ : i G N d } to be another i.i.d. collection of positive 
integer valued random variables with distribution function F(— 1 — G) and inde- 
pendent of {Xi : i G N d }. Let C := U{ ieN d| Xi =i}(i + [0, Pi] d )- We first consider 
eventual coverage of N d by X\ . 

Proposition 3.1. Let d > 2 and < p < 1. 

(a) if limine jG(j) = then P p (C eventually covers N d ) = 0, 

(b) i/liminfj^oo jG(j) > then P p (C eventually covers N d ) = 1. 

We sketch the proof for d = 2. For i,j G N let A(i,j) := & C}. Clearly, 

P(A(k,j)nA(i,j)) = P(A(k - i,j))P(A(i,j)) for k > i, 

i.e., for each fixed j the event A(i,j) is a renewal event. Thus, if, for every j > 1, 
P(A(i, j)) = oo then, on every line {y — j}, j > 1, we have infinitely many 
z's for which is uncovered with probability one and hence N rf can never be 

eventually covered. 

To calculate P p (A(i,j)) we divide the rectangle x as in Figure 1. For 
any point (fc, I), 1 < k < i — j and 1 < / < j, in the shaded region of Figure 1, we 
ensure that either Xtk,l) = or P(k,i) < k +j — 1. The remaining square region in 
Figure 1 is decomposed into j sub squares of length t, 1 < t < j — 1 and we ensure 
that for each point (fc, I) on the section of the boundary of the sub square t given 
by the dotted lines either Xr^n — or p(k,l) < t. So, 

j—l i—j 

P p (A(i,j)) = (l-p)l[(l-p + pF(t-l)) 2t+1 l[(l-p + pF(k + j-l)y 
t=l fe=l 

= (i- P )nV-pG(t)) 2t+i n (i-pg(w- (3) 

t=i fc=i+i 

Now choose e > such that pje < 1 and get N such that, for all i > N, iG(i) < e. 
Taking Cj := n?=J (1 ~ pG{t)) 2t+1 from © we have that 



Y^PviMhl)) = (l-p^Y^W^-pGik + j)) 3 

-N i=N k=l 

oo 

= { l -P)Cj^Z^ (say)- (4) 



i—N 
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J, 3) 



— • 



Fig 1. Division of the rectangle formed by [l,i] X [l,j]. 



For m > N we have 



(1 -pG(m + l)) j 



> 



1 - 



pe 



1 



pe 



»-p) 



fc=2 



kl 3 
k 



(m + l) fe 



1 - 



Pj' e , ff(ra,p,j, e) 
n + 1 (m + l) 2 



for some function g(m,p,j,e) bounded in m. Thus by Gauss' test, as pje < 1 we 
have Y^Ln ei = 00 an< ^ nence -fp(^(*> i)) — 00 • This completes the proof of 

the first part of the proposition. 

For the next part we fix r\ > such that 77 < lim infj_>oo jG(j) and get N\ 
such that for all i > Ni we have «G(i) > 77. Also, fix < p < 1 and choose a 
such that < exp(— prj) < a < 1. Let A/2 be such that for all j > N2 we have 
(1 - p-qj^ 1 )! < a. For N := max{iVi, N 2 }, let i,j e N be such that j > N and 
£ > j. Define := {(z, j) ^ C}. As in @ we have 



j'-i 



p p (A(i, j)) = a -p) n a - poo- + n (1 pgw) 2 



(5) 



fc=i 



t=i 
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Taking Cj := Yit=i (•"■ -pG 1 ^)) 2 '" 1 " 1 , we have from (0 and our choice of j, 

oo oo i—j 

J2P P (A(iJ)) = (l-picj^lKl-pGik + j)? 

i=N fe=i 

OO 



oo 



For m > N 



Jm+l 



i=N 

(l-pG(m + l)y 

< [1-P V J 



m + 1, 
pj'77 h(m,p,j,T)) 



m + 1 (m + l) 2 



(6) 



for some function h(m,p, j,r]) bounded in to; thus by Gauss' test, if > 1 then 

X^at &» < 00 an d hence Pp(Mhj)) < 00 . 

Now, for a given p, let j' := sup{j : pjrj < 1} and jo := max{j' + 1,N}. We 
next show that the region Qj := {(11,22) € N d : ii,Z2 > jo} has at most finitely 
many points that are not covered by C almost surely; there by proving that C 
eventually covers N d . For this we apply Borel-Cantelli lemma after showing that 
^(ii,«)efl Jo P P (Mh,i2)) < 00. 

Towards this end we have 

00 /io+tre-1 

=2(i- P )y( n (!-p G w) 2t+i 

00 k—rn \ 

x £ J](l-pG(j + rn + i))A+m 

fc=m+l t=l / 

00 jo + fe-1 



+e n (i-pow)^ 1 . 



fe=i t=i 



Observe that 



00 k — m 

@m •— Yi n o- - p g ^° +m+ i )) jo+m 

fc— m+1 «— 1 

hence as in © and the subsequent application of Gauss' test, we have that, for 
every to > 1, a m < 00. 

Now let 7 m := n^Li™" (1 — pG(t)) 2t+1 a m . Note that an application of the ratio 
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test yields Em=i 7™ < 00 > m deed from Q, 

7m+l 
7m 

_ (1 - pG(io + m)) 2j0+2m+1 EZi n:=i(l ~ pG(jo + m + 1 + i))*>+m+i 



Jo+m 



(1 - pG(io + m + l)) JO+m 1 + Es= 2 ~ PGO'o + ™ + 0) 

<(1 pG(j +m)) 1 + Er=i n-=i(l - PG(J0 + m + 1 + *))io + ™ " 

Since <r m < 00 for all m > 1, both the numerator and the denominator in the 
fraction above are finite. Moreover, each term in the sum of the numerator is less 
than the corresponding term in the sum of the denominator; yielding that the 
fraction is at most 1. Hence, for < a < 1 as chosen earlier 

< (l-pG(j +m)Y" +m+1 

7m 

< a. 

This shows that Em=i 7™ 00 ano - completes the proof of part (b) of the propo- 
sition. 

It may now be seen easily that p u and p l satisfy the conditions of Proposition 
3.1 and thus Theorem 1 1 . 41 holds . 



4. Markov Model 



The relation between the Poisson model and the discrete model explained in Section 
3 shows that Theorem l 1 . 31 would follow once we establish Theorem ll.5l In the setup 
of the Theorem Ol for each k 6 N let A k := {k g C}. To prove Theorem ILTT a) 
we show that ^2 k P(A k ) < 00 ano - an application of the Borel-Cantelli lemma 
yields the result, while to prove Theorem 11.51 (b) we show that ^ k P{A k ) = 00. 
However, the A k 's are not independent and hence Borel-Cantelli lemma cannot be 
applied. Nonetheless using the Markov property one can show that P(A k fl Ai) = 
P(Ak-i)P(Ai) and therefore, Aj's are renewal events; so by the renewal theorem, 
if EiSj P(Ai) — 00 then Ai occurs for infinitely many i's with probability one. 

For k > 1, let P (A k ) = P(A k \ X 1 = 0) and P 1 (A k ) = P(A k \ X 1 = 1). The 
following recurrence relations may be easily verified 

P (A k+1 ) = p 00 Po(A k ) + poiPi(A k ) (7) 
Pi(A k+1 ) = F(k - 1) [ Pw P a (A k ) +p n P 1 (A k )} . (8) 

We use this to prove Theorem ll.5f b) first. Let ^o( s ) = J2kL ko Po(A k )s k and 
^i( s ) = Efcfco Pi(A k )s k denote the generating functions of the sequences 
{Po{A k ) : k > ko} and {P\(A k ) : k > k } respectively, where ko is such that for 
a given e > and C — L + e > (where L is as in the statement of the theorem), 
fc +(l-C)>0, P (A kQ )>0, Px{A ko )>0, and F(fc-l) > for k > k . 

Such a ko exists by the condition of the theorem. 

Using the recurrence relations (JJJ and JSJ) we obtain 



*i(*)P(s) >Q(s)B(s) + R(s), 



(9) 
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where 



P(s) 


= (1 




= (1 


Q(s) 


= (1 


R(s) 


= (1 



s) 2 (l - C)pn + (1 - C)p w pois(l - p 00 s] 
+Piospoi(l - Poos) +PwPoaPois 2 

D feo-l TDl A, \ _1_ (U~ _1_ 1 _ ^ W .. c^O / 



+ Pw s ko+1 pooP (A ko ) 
From @ we have for any < t < 1 



Jo 



P(s) 



Now for s < 1, §W = h — h i — n — for some real numbers 

' P(s) 1-poos 1-s l-s(l-poi-pio) ' 

D, E, F. It may now be seen that ^1(1) = 00 whenever E > 0. Also the recurrence 
relations show that Wo(l) = 00 whenever Vt'i(l) = 00. A simple calculation now 
yields that E > if and only if — ^ — < i = Since e is arbitrary, we obtain 

J ■' pio+poi C L+ti •" 

Theorem IQTb). 

The proof of Thcorcm ll.Sf a) is similar. 
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